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Abstract. Via a Cole-Hopf transformation, the multicomponent linear heat hierarchy leads 
00 ■ to a multicomponent Burgers hierarchy. We show in particular that any solution of the latter 

also solves a corresponding multicomponent (potential) KP hierarchy. A generalization of 
the Cole-Hopf transformation leads to a more general relation between the multicomponent 
' linear heat hierarchy and the multicomponent KP hierarchy. From this results a construction 

r*j ■ of exact solutions of the latter via a matrix linear system. 
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1 Introduction 

The well-known Cole-Hopf transformation <fi = il) x ih~ l translates the nonlinear Burgers equation 
into the linear heat equation [H [2J [3] . This extends to a corresponding relation between the 
Burgers hierarchy and the linear heat hierarchy, and moreover to their matrix generalizations. 
The Cole-Hopf transformation also generates solutions of the KP hierarchy from (invertible) 
solutions of the linear heat hierarchy [H (U [6], and this extends to the corresponding matrix 
hierarchies [3 El [9] . There is also a generalization of the Cole-Hopf transformation that produces 
a solution of the (scalar or matrix) KP hierarchy from two solutions of the matrix linear heat 
hierarchy, connected by an additional relation [7J[8]. Furthermore, this includes a construction 
of (matrix) KP solutions from a matrix linear system (see also 111]), which may be regarded 
as a finite-dimensional version of the Sato theory [121 113j . In this work we extend these results 
to the multicomponent case. Admittedly, this is not a very difficult task, on the basis of our 
previous results [Tl]. We also take this opportunity, however, to present these results in 
a concise form and to add some relevant remarks and examples. 

The usual multicomponent KP (mcKP) hierarchy contains subhierarchies that generalize 
the matrix KP hierarchy by modifying it with some constant matrix B different from the 
identity matrix. We are particularly interested in such matrix hierarchies different from the 
ordinary matrix KP hierarchy. The Davey-Stewartson equation |14t 115] , a two-dimensional 
nonlinear Schrodinger equation that appeared as a shallow-water limit of the Benney-Roskes 
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equation [T7] , emerges from such a modified matrix hierarchy (see Section [3|) . It is of course 
well-known to arise from the two-component KP hierarchy. 

Section [2] recalls the Cole-Hopf transformation for a multicomponent Burgers (mcBurgers) 
hierarchy and Section [3] reveals relations between the mcBurgers and the mcKP hierarchy. Sec- 
tion U] presents the abovementioned kind of generalization of the Cole-Hopf transformation which 
determines solutions of the mcKP hierarchy. From this derives a fairly simple method to con- 
struct mcKP solutions. This is the subject of Section which also presents some examples. 
All this is closely related to a multicomponent version of a matrix Riccati hierarchy, as briefly 
explained in Section Section [7] contains some conclusions. 

In the following, let A be an associative (and typically noncommutative) algebra over the 
complex (or real) numbers, with identity element I and supplied with a structure that allows 
to define derivatives with respect to real variables. Let B be a finite set of mutually commuting 
elements of A. Although this assumption is sufficient to establish the results in this work, further 
assumptions should be placed on B in particular in order to diminish redundancy, see the remark 
at the end of Section [2j With each B G B we associate a sequence of (real or complex) variables 
t_B = (£b,i, • • • -,tB,ni ■ ■ •)• Furthermore, for a function / of {tsj-sgg we define a Miwa shift with 
respect to B £ B by /[A] s (- • • ,ts, . . .) = /(..., t B + [A],...), where [A] = (A, A/2, A/3, . . .) and A 
is an indeterminate. ft n denotes the partial derivative of / with respect to the variable t n . 



2 Cole-Hopf transformation 

for a multicomponent Burgers hierarchy 

Let us consider the multicomponent linear heat hierarchy 

lH B ,n = B n d n (iP) VBeB, n = l,2,..., (2.1) 

where d = d x is the operator of partial differentiation with respect to a variable x, and ip has 
values in .4. Since ipt B i = B^ x , this implies the ordinary linear heat hierarchy ipt B „ — &t B l (VO- 
Any two flows (|2.ip commute as a consequence of our assumptions for B (the elements do not 
depend on {ts}BeB an d commute with each other). A functional representation^ of (12. ip is 
given by 

\- l ^-^ [x]B ) = B^ x VBeB. (2.2) 

Proposition 1. If ip is an invertible solution of the above multicomponent linear heat hierarchy, 
then 

<t> = (2.3) 

solves the multicomponent Burgers (mcBurgers) hierarchy associated with B, given by the func- 
tional representation 

n B (<P,\) = VBeB, (2.4) 

where 

fisOk A) := A" 1 ^ - 0_ [A]B ) - (B0 - </>. [x]B B)<f> - B^ x . (2.5) 



1 We use this (not quite satisfactory) term for an equation that generates a sequence of equations by expansion 
in powers of an indeterminate. 
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Proof. We have to consider the following system, 



The integrabihty condition (ip x )-[\] B = (V'-[a]b)z yields (|2.4p . The further integrabihty condi- 
tion (^_ [a]si )_ Mb2 = (V_ Ms2 )-[A] Sl is 



^(A- 1 ^ - 0-[A] Bl ) + 0-[A] Bl ^i0) = Bi^" 1 ^ - 0-M fla ) + ^-Mb 2 ^20), 



which is satisfied as a consequence of (j2.4j) and B2] =0. I 



(|2.3p is a Cole-Hopf transformation. The first equation that results from (|2.4p is 



w - 



which has been called C-integrable N-wave equation in [9]. If B = /, this reduces to ij x = x. 
But it is a nontrivial nonlinear equation if </> does not commute with B. The next equation that 
results from (|2.4p is the (noncommutative) Burgers equation 

Remark 1. In order to avoid redundancy and to maximally extend the hierarchies, further 
conditions have to be imposed on B. In particular, the elements of B should be linearly indepen- 
dent, since linear combinations correspond to linear combinations of hierarchy equations. But 
since products of elements of B also generate (independent or redundant) commuting flows, the 
problem is more subtle. If the algebra is semisimple, it admits a maximal set of commuting mu- 
tually annihilating idempotents E a , a = 1,...,N, hence E a Ei = 8 a ^E a for all a, b = 1, . . . , N. 
Then B = {E a }% = 1 is an optimal choice. In fact, in the following we do not really address those 
cases of (non-semisimple) algebras where such a choice does not exist. Rather it turns out that 
some more flexibility in the choice of B can be used to obtain certain integrable systems within 
this framework in a more direct way, see Section [3j 

3 Multicomponent KP and relations with the multicomponent 
Burgers hierarchy 

For B e B let 

£ B {\) :=I-X{uj B {X) + Bd). 
The "discrete" zero curvature condition 



(3.1) 



then leads to the two equations 



A (uBaifi) -u B ^)-[\] Bl ) +wb 2 (/x)_ [a]bi o;b 1 (A) + B2U Bl (\) x 

= M~Vi?i(A) -UB 1 (X)-\ M ] B )+WB 1 W-\ M ] B WB a (p)+B 1 V Ba (p) 



x 



and 



B 2 u} Bl (\) - u Bi (X)_ Mb2 B 2 = BxWBzili) - wb 2 (m)-[a] Si £i- 



The last equation is solved by 



u B (X) = Bcf > -^ [x]B B, 
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and the first equation can then be written in terms of (|2.5p as 

B 2 n Bl (0, a) - n Bl o, \)- Mb b 2 = B 1 n B2 (0, n) - n B2 (<P, aO-m^i vb 1} b 2 <=b. (3.2) 

We take this as our defining equations of the (more precisely, potential) multicomponent KP 
(mcKP) hierarchy associated with £>H 

Remark 2. Choosing B\ = B 2 = B in (|3.2p and summing the resulting equation three times 
with cyclically permuted indeterminates, leads to 

3 

(KHt - 0-Mb) + ^ { x t ] B Bcb)_ [Xj]B 5 = 0, 

i,j,k=l 

which is a special case of the functional form of the mcKP hierarchy in [27| . 

Let us take a closer look at (|3.2p with B\ = B 2 = B. Its A-independent part is 

Bn B {<t>, o) - n B {<t>, o)_ Mb b = [b, n B (<p, M )], 

where £l B ((j), 0) = </>< s x — [B, (f>\<j) — B(j) x . To first order in (i this gives 

Bfo B xX B - l{B,<j> tB ltB1 ] + \[B,4> tBa ] = B(j>t B1 [B,(j)] - [B, 4>]4>t B1 B. (3.3) 

(|3.2p is the integrability condition of 

n B (0,A) = S0-fl_ [A]B B VBeB, (3.4) 

with a new dependent variable 9. (|3.4|) represents the mcKP hierarchy in terms of two dependent 
variables. By comparison with (|2,4p . this has the form of an inhomogeneous mcBurgers hierarchy. 
The following is an immediate consequenc^l. 

Proposition 2. Any solution of the mcBurgers hierarchy also solves the mcKP hierarchy. 

Proof. ([33D becomes {23]) if 9 = 0. ■ 

The representation (|3.4p of the mcKP hierarchy has the advantage that each equation only 
involves a single element from the set B. To order A , (|3.4|) yields 



[B, </>]<!>= [B,9]. (3.5) 



If B = I, then (|3.5|) reduces to <^ s t = (f> x - Otherwise this is a non-trivial nonlinear equation. 
Subtracting (|3.5p from (|3,4p . leads to 

A- 1 ^ - 0_ [A]fl ) - iSjl - (0 - 0_r A]fl W = (0 - 0-[a]b) b - (3-6) 

The two equations (|3.5|) and (|3.6p are equivalent to (|3.4|) . (|3.6|) does not involve derivatives with 
respect to x. For fixed B £ B, it represents the KP hierarchy in A [7J, with product modified 
by B. 



2 See also e.g. [T8l[Tl[20l[2Tl[23|l[23[M[25l[26l[2^ different formulations of such a multicomponent KP 
hierarchy. We should also mention that the conditions imposed on the set B can be relaxed while keeping the 
hierarchy property, see [26] for example. 

3 It was first noted in [4] that any solution of the first two equations of the (scalar) Burgers hierarchy also solves 
the (scalar potential) KP equation. In [281 129] the (first two) Burgers hierarchy equations have been recovered 
via a symmetry constraint from the KP hierarchy and its linear system. 
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Proposition 3. For any B 6 B, as a consequence of the mcKP hierarchy (|3.4p . Bcj) and also 
4>B solves the ordinary (noncommutative) KP hierarchy. 



Proof. This is an immediate consequence of (|3.6p . ■ 
To first order in A -1 , (|3.6|) yields 

Differentiating (|3.5j) with respect to tB,x, and multiplying it by B from the right, we can use 
the last equation to eliminate 6. In this way we recover (13. 3f) . As shown in Example [2] below, 
(|3.3p generalizes the Davey-Stewartson (DS) equation [TBI CEH 13 Q2] . Eliminating 8 from (|3.6|) 
by use of (13, 5p thus leads to a (generalized) DS hierarchy. 

Example 1. Let B = a, where a 2 = I. Decomposing as 

<P = D + U, where D := |(0 + afta), U := - after), 

(|3.3p splits into the two equations^] 

D tl a-D x -2U 2 = 0, (3.7) 
U t2 a + U xtl +2{U,D tl } = 0, (3.8) 

where we write t\, ti instead of an d t_B,2- Let now A be the algebra of 2 x 2 matrices over C, 
and a = diag(l, — 1). Then D and U are diagonal and off-diagonal parts of <f>, respectively. 
Writing 

U=( ° " J, hence u:=0i 2 , u := 02,1, (3.9) 



* u , 
and introducing 

s := tr(0) = 0i 5 i + 02,2, r := 0i,i - 02,2, 
we obtain the system 

u t2 - u xtl - 2us tl = 0, v t2 + v xtl + 2vs tl = 0, (3.10) 

and sjj = r x , s x = r tl — Auv. The integr ability conditions of the latter two equations are 

st!t! ~ s xx = 4(uv) x , r tltl - r xx = 4(uv) tl . (3.11) 

The two equations (I3.10p together with the first of (13. lip constitute a fairly simple system 
of three nonlinear coupled equations, where all variables can be taken to be real. Allowing 
complex dependent and independent variables, after a complex transformation the system for 
the dependent variables u, v, s can be further reduced to the DS equation, see the next example. 
Of course, the above system (]3.10l) and (|3.1ip can also be derived from the usual two-component 
KP hierarchy (see e.g. [27J), and the transformation to DS is well-known. 

Setting v = 0, we obtain from (|3.10p and (|3.1ip the following linear equation^!, 

ut 2 ~ Uxh ~ 2l P u = 0, (ftiti ~ fxx = 0, (3.12) 

where <p := st L = tr(0)^ 1 . This is probably the simplest system that possesses dromion solutions, 
as observed in Example [3] in Section [5j 



4 A constant (with respect to ti) of integration has been set to zero in order to obtain (|3.7|) . The latter equation 
can be obtained more directly as the diagonal part of (|3,5l l. 

5 If u does not depend on i2, then the first equation is part of a Lax pair for the Nizhnik-Novikov-Veselov 
equation. 
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Example 2. We continue with the previous example and perform a transformation to the DS 
equation by first allowing the dependent variables to live in a noncommutative algebra. In this 
way we obtain a certain noncommutative generalization of the DS equation (see also |30l [3Tj for 
matrix DS versions). In terms of 

F = D tl + (3U 2 , 

with G C, HMD becomes 

F tl a -F x = -(3(U 2 ) X + [U 2 ) tl {2 + (5a). 

Differentiating this with respect to x and with respect to t\, respectively, and eliminating mixed 
derivatives of F from the resulting two equations, yields 

F tlh ~ F xx = -(3(U 2 ) XX + 2(U 2 ) xtl + (U 2 ) tltl (2a + (3). (3.13) 

Furthermore, fl3.8f) takes the form 

U t2 a + U xtl = -2{U, F} + A(3U 3 . (3.14) 

Let now A be the algebra of 2 x 2 matrices over some unital associative algebra with unit denoted 
by 1, and a = diag(l, —1). Using (|3.9|) and writing 

F = diag(/,fif), 

(|3.13|) and (|3.14|) result in the following equations, 

u t2 - u xtl = 2(fu + ug) - A/3uvu, 

v t2 + v xtl = -2(gv + vf) + A(3vuv, 

/titi - fxx = -(3(uv) xx + 2(uv) xtl + (2 + 0)(uv) tltl , 

dtxtx - 9xx = -P{vu) xx + 2(vu) xtl - (2 - f3)(vu)t ltl . 

In terms of the new variables 

1 + i 1 — i 

V= — ^ {x + ih), z = — -j=-[x - iii), (3.15) 

with a constant and with the choice (3 = i , this becomes 

ut 2 + u zz + —u yy = 2(fu + ug) - 4i uvu, 

-v t2 + v zz + \v yy = 2(gv + vf) + 4ivuv, 

. . . 1 . . 2i . , 
l + 2i ){uv) zz — 2^ uv >yy ( uv )yz, 

-1 + 2\)(VU) ZZ + -^(Vu)yy + —(Vu) yz . 

If the dependent variables take their values in a commutative algebra, then we obtain 

1 2 
ut 2 + u zz H — 2 u yy = ^P u ~ ^ u v ' 

1 2 

-v t2 + v zz + — v yy = 2pv + AlV u, 

K 

Pzz - -2Pyy = 4i (™) 22 , 
K 







fzz 


If 




1 


9zz - 





(3.16) 
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where 

p ■- / + g = ip + 2iuv (3.17) 
(with ip = tr(</>) tl ). Setting ^2 = — i * and 

v = i eu with e = ±1 
and the complex conjugate u of u, one recovers the DS system 

iu t + u zz + \u yy = 2pu + 4e\u\ 2 u, p zz - \p yy = -4e(\u\ 2 ) zz , (3.18) 

where now y, z, t, k 2 , p are taken to be real. Clearly, this system is more quickly obtained 
from (|3.10p and (|3.1ip by application of the transformation of independent variables given 
by (|3.15p and t<i = —it, and the transformation (|3.17p of dependent variables. But on our way 
we obtained the system (|3.16|) which (with £2 = —it) may be of interest as a noncommutative 
version of the DS equation (see [301 [31] f° r alternatives). 

We conclude that the DS equation (and a corresponding hierarchy) is obtained from a matrix 
KP hierarchy, but the latter has to be generalized by introduction of a matrix B different from 
the unit matrix. Of course, this hierarchy is embedded in the usual two-component KP hierarchy. 

3.1 The associated Sato— Wilson system and its translation 
into a Burgers hierarchy 



(|3.ip is the integrability condition of the linear system 

^_ [A]s =8 B {\)^ VBeB. 
If W is an invertible solution of 

£ B (\)W = W_ [x]b £ ,bW, So,bW : = 1 ~ XBd ' ( 3 - 19 ) 

then the linear system is mapped to ip-\\] g = £o,b(^)iP where ip := W~ x ip. The latter is the 
linear heat hierarchy (|2.2|) . The ansatz 

W = I + ^w n d~ n , 

n>\ 

inserted into the functional form (|3.19p of the Sato-Wilson equations, leads to 

X^iwn - w n ,-[a] s ) - Bw UtX - {B4> - 4>_\x\ B B)w n = Bw n+ i - w n+1 „ [x]B B, (3.20) 

where n = 0, 1, . . . and w$ = I. From the n = equation we get w\ = —(f). In terms of 9 := —W2, 
the n = 1 equation turns out to be the functional form (|3.4p of the mcKP hierarchy. Let us 
introduce 

$ := A T + ei((f), —W2, —W3, . . .), (3.21) 
where e\ = (J, 0, . . .), and 



A 



/ / 
0/0 



(3.22) 
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is the shift operator matrix with transpose A T . Then (|3.20p can be expressed as the following 
mcBurgers hierarchy, 

SIb(*,\) = \- 1 ($-$-[\] b )-(B$-§_ [x]b B)$-B$ x = VBeB. (3.23) 

The linear system of the mcKP hierarchy has thus been reformulated as an mcBurgers hierarchy 
with an infinite matrix variable of a special form. 

We note that (|3.2ip has the form of a companion matrix, a block of a Frobenius normal 
form matrix. This makes contact with recent work in [S], where several integrable systems were 
recovered from equations of a Burgers hierarchy with dependent variable of Frobenius normal 
form. All these systems are known to arise as reductions of the mcKP hierarchy. 

As a consequence of (|3.23j) with $ of the form (|3.2ip . it follows that <f> = ej$ei solves the 
mcKP hierarchy in A. 

In the following section, we shall see that $ of the form (|3.2ip results, as a particular case, 
from a quite general result about solutions of a somewhat generalized mcKP hierarchy. 

4 Generalization of the Cole— Hopf transformation 

We generalize (13. 4j) to 

n BtQ (<j>,X)=BQ-@_ [x]B B VBeB, (4.1) 

where 

B>Q ($, A) := A- x ($ - d>_ [A]B ) - (B$ - <D„ [a]bJ B)Q$ - B<S> X , (4.2) 

with a constant object Q. We assume that the constituents are elements of linear spaces and that 
the products are defined. (|4.ip determines a generalization of the mcKP hierarchy, since non- 
linear terms now involve Q (which modifies the product). The following generalizes a theorem 
in [7] (see also [8]). 

Theorem 1. Let X, Y be solutions of the multicomponent linear heat hierarchy, i.e. 

X~\X - X_ [X]B ) = BX X , \-\Y - Y_ [X]B ) = BY X , (4.3) 
for all B £ B, and furthermore 

X X = RX + QY, (4.4) 
with constant objects R, Q. If X is invertible and if all B G B commute with R, then 

$ = YX- 1 (4.5) 

solves the mcKPQ hierarchy (|4.1|) with = QR. 

Proof. Using (|4.5p in the definition (I4.2p . we have 

n B) g($, A) = [B$> - $ ax]b B)(X x - QY)X~ 1 + (X'^Y - Y_ [a]b ) - BY X )X' X 
~ §- { x ]B (X\X - X_ [X]B ) - BX X )X~\ 

which reduces to 

n B , Q (^,X) = (B^-^_ [x]B B)R 
as a consequence of ([4T3]) and (fO|) . Since [B, R] = 0, this takes the form (|4TT|) with 6 = $R. ■ 
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Now we set up the stage for applications of the theorem. Let A(M, N) := Mat(M x N, C)®A, 
where Mat(M x N, C) is the space of complex M x N matrices. Let 0, Y take values in 
A(M,N), and X in A(N,N). Furthermore, let Q G A(N,M) and R G A(N,N) commute with 
all BeB. 

If Q has ranfc one over A, in the sense that Q = VU J with constant vectors U, V, with entries 
in A, and if J7 and V commute with all B £ B, then = £/ T< l ) V solves the mcKP hierarchy in A, 
provided that <3? solves (|4.ip . In this way, any solution X, Y of the linear equations formulated 
in the above theorem generates an A- valued solution <j) of the mcKP hierarchy (13. 2\\ . 

Choosing M = N, (|4,5p is a Cole-Hopf transformation if Y = X x . Then (|4.4p becomes 
(/ — Q) X x = RX. Let efe denote the X-component vector with all entries zero except for the 
identity element in the kth row. 

N-l 

a) Setting Q = e^ej^ and R = e fc e fc+i (which is the left shift operator: Re^ = e^-i, 

fc=i 

k = 2, . . . , N, and i?ei = 0), one finds that (|4.4p restricts X to the form of a Wronski matrix 
(see also [8]), 



d(xW) 



d(xM) 



d(xW) 



are independent functions with values in ^4. The remaining assumption 
in the theorem requires them to be solutions of the multicomponent heat hierarchy. 

b) Let R = A T with the infinite shift operator matrix (|3.22p . and Q = e\ej. Then (|4.4p says 
that X has to be a pseudo- Wronski matrix 



X 



( X« 

V 



X( 2 ) 

9 -2 x (2) a -2^(3) 



\ 



(4.6) 



where d~ l is the formal inverse of d. This structure appeared in [9] (equations (14) and (50) 
therein). (|4.3p demands that X^ 1 ) , X^ 2 ) , . . . solve the multicomponent heat hierarchy. With (|4.6p . 
$ = X^X" 1 has the form (|3.2ip and hence determines a solution of (|3.23p . 



5 Solutions of the multicomponent KP hierarchy 
from a matrix linear system 

In order to derive some classes of mcKP solutions via theorem [1] more explicitly, in the framework 
specified in Section [4] (after the theorem) we extend (|4.4|) to 



Z x — HZ, 



(5.1) 



where 



Z 



X 
Y 



H 



R Q 
S L 



with new constant objects L G A(M, M) and 5 G A(M, X) that commute with all B G B. The 
two equations (14, 3h then combine to 



\-'(Z - Z_ [X]B ) = BZ X V B e B. 
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Taking (|5.ip into account, this is equivalent to 

Z tBn = {BH) n Z V -B G £>, n = l,2,.... (5.2) 

Note that B and H commute as a consequence of our assumptions. With a suitable choice of 
the algebra A, the solution of the matrix linear system (|5.ip . (|5.2p is given by 

Z(x,t) =e^ x ' t ' H ' B ^Z , (5.3) 

where 

oo 

t; , fl) := xtf + J] £ t B , n (HB) n , 
BeBn=l 

and t stands for {t^}_BeB- Decomposing Z into X and Y, theorem [T] implies that = 
solves the hicKPq hierarchy (|4.ip . Furthermore, if rank(Q) = 1 over A, hence Q = VU T with 
constant vectors U and V, then the A- valued variable 

(f, = U J $V 

solves the corresponding mcKP hierarchy. 

The exponential in (|5.3p can be computed explicitly if additional assumptions are made 
concerning the form of H (see |llj . in particular). Then is obtained via (|4.5p . 

Case 1. Let S = and 

Q = RK - KL (5.4) 
with a constant N x M matrix K (over .A) that commutes with all B G B. Then we obtain 

$ = e5 (,,t;L,B) Ce - f (x,t;R,iJ)( /jv _ ^(^t^B^-^t^B)) (g 

where /jv is the N x N unit matrix over A (so that the diagonal entries are the identity I 
in A), and C is an arbitrary constant M x N matrix (with entries in A). $ solves the ulcKPq 
hierarchy (associated with B), with Q given by (|5.4p Fl If moreover Q = VU T with vectors U, V 
that commute with all B G B, then <p = U T $V solves the mcKP hierarchy in A. Of course, it 
remains to solve the rank one condition (over A) 

RK - KL = VU T . (5.6) 

If M = N and if C is invertible, then (15.5f) simplifies to 

$ _ f e $(x,t-,R,B) c -l e -$(x,t-,L,B) _ K)~ X , 

which remains a solution if we replace C _1 by an arbitrary constant N x N matrix C (with 
entries in A). 

Example 3. Choosing the components of the matrices L, R as 

Lij = PiSijI, Rki = qkfall, (5-7) 



In particular, if M = N and Q = In, then $ (with K, L, R solving RK — KL = In) is a solution of the 
N x N matrix (over „4) mcKP hierarchy. 
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with constants pi, q^, (|5.6p is solved by 

K k j = — - — u k Vj, (5.8) 
Qk ~ Pj 

where u& and Vj are the components of U and V, respectively. We elaborate one of the simplest 
cases in some detail. Let us choose A as the algebra of 2 x 2 matrices over C, L = pfy, R = ql2, 
with constants p and q, Q = I = I2, and B = {B} with B = diag(l, —1) (motivated by the 
examples in Section [3]). Then we have £(x,t\ L, B) = diag(£+(a;,t;p),£_(x,t;p)) with 

00 00 

e±(x,t;p) :=px + Y,P 2n t2n±Y,P 2n+H ^+l ^ 

n=l n=0 

(writing t n instead of tB,n), and £(a:,t; R, B) is obtained by exchanging p by q in these expres- 
sions. Writing 

C=( Cl C2 ) (5.10) 
with constants a, we obtain (with U = V = i^jj 

1 / Cie t+(x,t;p)-Z+(x,t;q) +f( X ,t) C2e e+(x,t;p)-C_(x,t;g) \ 



where 



._ C1C4 - C 2 C 3 f + ( a!l t3»)+f-(*,ta»)-e + (!r,t;g)-€-(a!,t; g ) 



/(x,t) — e 

p-q 

V(x,t) -.= 1 + - J—fdgS+Ca.tipJ-C+fo**) + C4e «-(^t;p)-?-(^;'?) + /( x t )). 

This is a solution of (|3.3|) . with .B = diag(l, — 1), and its hierarchy, and its components thus 
provide us with the solution 

u = C Jl e «+(M;p)-C-(z,t;g) v = — gf- 

s = i( Cie €+(*,t»)-e+(a!,tiff) +C4 gC-(x,t»»)-e-(*,t;«) +2 /) (5.12) 

of the system (|3.1U|) . (|3.11|) . is regular (for all t) in particular if all constants are real, 
C1C4 > C2C3, and either p > ci > 0, C4 > 0, or p < q, c\ < 0, C4 < 0. Fig. [1] presents 
a dromioio solution within this family. 

If C3 = 0, we have v = and (|5.12p determines a solution of the linear equations (|3.12p l^l An 
extremum of u for a regular solution then moves (in "time" £2) with constant amplitude along 
the curve given by 

x = -(p + q )t 2 -—— — -log ( ClC4 2 J , h = — ^log ( --J— J ■ 

2{p-q) \{p-q) J p-q V q\ c iJ 

The last expression shows that, for a dromion solution, p and q must have opposite signs. 



7 Since £±(x,t;p) — £±(x,t;q) = (p — q)(x ± ti) + (p 2 — q 2 )t2 + • ■ • , this solution becomes independent of £2 
(i.e. "static") if p = — q. 

8 The characteristic features of a dromion are its exponential localization and that it is accompanied by a field 
structure of intersecting line solitons. See [321 [33] 1341 1351 1361 1371 1381 1391 1401 141| for the Davey-Stewartson case, 
and especially [32] f° r an illuminating structural analysis and the appearance of dromions as solutions of other 
equations. 

9 In this case, T> factorizes, 

X>= [1+ Cl gM^irt-M^t^A [ 1 4. 04 e «-(x,t;p)-{_(x,t;g) 

p-q J \ p-q 
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Figure 1. A dromion solution of (|3 . 10[) and (|3 . 11 [> at *2 = 0, given by (|5.12[) with p = 2, q = — 1, 
c\ = C4 = 1, C2 = 1/2 and C3 = 0. The left plot shows <p = s tl — tr(0) tl . As a consequence of 
C3 = 0, we have v = 0. Hence this is actually a solution of the linear equations (|3.12[) and ip solves the 
linear wave equation. Regarding t\ as an evolution parameter, the plot of (p shows two colliding humps 
(with amplitudes having opposite signs) that annihilate at t\ = 0. With 7^ C3 < 2, the plots remain 
qualitatively the same as long as C3 is sufficiently far below the upper bound, and v attains a shape 
similar to that of u. 




Figure 2. Solutions of (j31Q)) and (f3TT]) at t 2 = 0, given by ([5TT]) with p = 2, a = c 4 = 1, c 2 = c 3 = 1/2. 
The first two plots, where q = 0, show kinks (see (|5.13[) V In the last two plots, where q = —1/5, these 
become exponentially localized structures. 



For q = in (|5.12|) . setting t n = for n > 2, we have 

u = c 2 p 2 {c lP + c 4 pe- 2ph + e - ptl (p 2 e- p( - x+pt ^ + (cic 4 - c 2 c 3 )e p{x+pt2) ) ~\ 

v = c 3P 2 (c iP + c lP e 2ptl + e ptl (p 2 e- p{x+pt2) + (cic 4 - c 2 c 3 )e p{x+pt2) ) ~\ (5.13) 

If C1C4 > C2C3 and c\p > 0, C4P > (or C1C4 < c 2 c 3 and cip < 0, c 4 p < 0), these functions 
obviously represent wedge-shaped kinks in the xii-plane, see also Fig. [2j Choosing p > and 
switching on a negative q, these wedges become localized and, around certain negative values 
of q, then take the dromion form. 

Fig. [3] shows plots of a two-dromion solution determined by (|5.7p with L = diag(3/2j 2/2)) 
R = diag(-2/ 2 , -(3/2)/ 2 ), anc0 

/ 1 i \ 
10 
~ 2 1 
V 5 3 / 

The diagonal 2x2 blocks of these matrices correspond to matrix data of single dromions. 
Such a superposition is obtained for any two solutions, provided that off-diagonal blocks of the 

10 If all lower-diagonal entries of C are zero, we obtain v — and thus a solution of the linear equations (|3.12[) . 
The plots are surprisingly insensitive with respect to changes in this range of parameters, as long as the off- 
diagonal entries in a diagonal block are not all close to zero and the determinant of the block is not close to 
zero. 
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Figure 3. A two-dromion solution of (|3.10|) and (|3 . at t?, = 1, determined by the data specified in 
Example [3j The left plot shows <p = = tr(</>) tl . 




Figure 4. A dromion solution of the Davey-Stewartson-I equation with e = — 1 at t-x = 0, see the end of 
Example [3] Here we chose pi = p 2 = (1 +i)/V2, qi = — <?2 = (1 — i )/V2, c\ = — c 2 = i \/2, C3 = C4 = s/2. 
Here p is given by (|3. lTf) . 

matrix K exist such that (|5.6p can be satisfied. This is so in the restricted case considered above 
(see (|5.7p and (|5.8p ). which in particular leads to multi-dromion solutions. Introducing non-zero 
constants in the off-diagonal blocks of C, leads to solutions with more complicated behaviour. 

Setting t n = 0, n > 2, the transition to the Davey-Stewartson system (13.18P with k = 1, 
which is the DS-I case, is given by the transformation of independent variables 

x = -^^(z-iy), h = -^-^(z + iy), t 2 = -it. 

The dependent variables are u and p, the latter given by (|3.17p . We have to take the additional 
constraint v = ieu into account (see Example [2]). One recovers a DS-I dromion within the class 
of solutions restricted by q± = p±, q 2 = —p2, c\ imaginary and C4 real, and C3 = ±i c 2 with sign 
corresponding to e = ±1. Fig. [5] shows an example. 

Case 2. Let M = N, R = L, S = 0, and 

Q = J+[L,K], (5.14) 

with constant N x N matrices K and J (over .A) that commute with all B S B. Furthermore, 
J has to commute with L, i.e. [J, L] = 0. Then 

$ = e Z(x,t;L,B) Ce -£(x,t;L,B) ^ + t; LB y_ K ^{x,t;L^) Ce -({x,t;L^ ^ ^ 

where C is an arbitrary constant N x N matrix (with entries in .4) and 

e(x, t;L,B):=x+^2Yl nt B , n L n ~ l B n . 
BeBn=l 
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If also Q = VU J with vectors U, V that commute with all B 6 B, then (j) = U J <&V solves the 
mcKP hierarchy in A. It remains to solve 

J + [L, K] = VU T . (5.16) 

A natural choice for J is the unit matrix In, but there are others. (|5.15p can also be written as 

$= ( f £(*,t;L,B)c e -t(x,t;L,B) + f'( Xj t; L, B) J - K)' 1 (5.17) 

with an arbitrary constant N x N matrix C@ If C is chosen such that it commutes with L 
and B, then $ and the corresponding solution <j) of the mcKP hierarchy are purely rational 
functions of the independent variables. A localized solution of this kind, hence with rational 
decay, is usually called a "lump". The following example in particular demonstrates that there 
can be weaker conditions that lead to solutions with rational decay. 

Example 4. Choosing B = diag(l,— 1), L = diag(pi,j>2) and Q = I 2 , (|5.14p is solved by 
K = diag(fei, k 2 ). Expressing C again as in (|5.10p . we find 



ci c 2 



e £+0z,t;pi)-§_O,t;p 2 ) 



C3e -f+0s>t;pi)+€-0B>t;P2) C4 



e e(x,t;L ) B) c , e -e(x,t;i ) B) _ 

with £±(x,t;p) defined in (|5.9p . and £'(x,t;L,B) = diag(^' + (x,t;pi),S i '_(x,t;p 2 )), where 

00 00 
£±(x,t;p) := x + J2 2nt 2 nP 2n ~ l ± J^(2n + l)i 2 n+ip 2 



n=l n=0 



Then (15.150 leads to the following solution of (13. 3D . with B = diag(l, —1), and its hierarchy, 
1 ( c l + d{£'_(x,t;p2)-k 2 ) C2e M^;pi)-^-ix,t; P2 ) 



^ X,t) V( X ,t) V C3e -5 + (-,t; Pl )+ 5 -(x,t; P2 ) C4 + d(eV(^t;pi)-Jfei; 

where d := C1C4 — C2C3 and 

£>(x,t) := 1 + c 1 (£' + (x,t;p 1 ) - h) + c 4 (C(x,t;p 2 ) - k 2 ) 
+ d(£' + (x,t;pi) - k 1 )(C(x,t;p 2 ) - k 2 ). 

The components u = 4>\ )2 and v = (f> 2 ,i, together with 

s = tr(0) = - (s, t;pi) + ^ (s, t; p 2 ) + c) , 

with a constant c, thus solve the system ()3.10p . (|3.1ip . For C3 = 0, this determines a solution of 
the linear equations (|3.12p . 

The transition to the Davey-Stewartson system (|3.18p with k = i , which is the DS-II case, 
involves the transformation of independent variables 

% = \^j={y + z), tl = \^ <yV ~ Z ^ h = -it. 

We set t n = 0, n > 2, in the following. The lump solution of the DS-II equation [43 1 144^ [TT] (see 
also |45 l l46 1l4"0ll47| ) is obtained as follows from the above formula. Besides taking account of the 
constraint v = ieu, we have to arrange in particular that the exponential in u becomes a phase 



11 In the transition from (|5.15p to (|5.17[) . one assumes that C is invertible with inverse C. But C need not be 
invertible in order that (|5.17|l determines a solution of the mcKPQ hierarchy. 
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factor (up to some constant factor), i.e. the real part of its exponent has to be constant. This 
requires setting p 2 = — ijPi- It is more tricky to find conditions on the remaining parameters 
such that T> ^ for all y, z, t, so that the solution is regular. Choosing 

c 4 = -ici, c 3 = ic 2 , k 2 = -iki, 

and renaming k±, p\ to k, p, we find that 

l/i |2 , ni |2\// ,\2 , / ,\2\ , 2 1 C2 1 2 



V = -(|ci| 2 + 2\c 2 \ 2 ){(y -y - v y t) 2 + (z - z - v z tf) + 



where = — 2Im(p), t> z = 2Re(p), and 

2Re(ci) 2Im(ci) 



yo = Re(fc) - | — „ — z = Im(fc) + 



|C1| 2 + 2| C2 |2' ~ U ""^ | Cl |2 + 2| C2 |2- 

The resulting DS-II solution 

u= C 2l e i(lm(p)y+R C (p)z-Rc(p 2 )t) 

V 

p = ^2 (d c il 2 + 2 N 2 ) 2 ((« - z o - M) 2 -(y-yo- v y t) 2 ) - 8\c 2 \ 2 ), 

with p defined in (|3.17p and e = 1, is regular whenever c 2 ^ and reproduces a well-known 
lump solution. 

Again, (lump) solutions can be superposed by taking matrix data of (lump) solutions as 
diagonal blocks of larger matrices L and C. It then essentially remains to determine the off- 
diagonal blocks of the new matrix K so that (15. 16ft holds. 



6 The matrix Riccati system associated 
with the multicomponent KP hierarchy 

Writing 

j^ri _. f Rn Qn \ 
\ S n L n J 

the matrix linear system (15.11) . (15. 2p implies the matrix Riccati system 

<$> x = S + L$-<f>R-$Q<f>, (6.1) 
$ tB n = B n S n + B n L n $ - $>B n R n - $B n Q n <Z>. (6.2) 

The two solution families presented in Section [5] solve this matrix Riccati system, with the 
respective conditions imposed on the matrices L, Q, R, S. 

Abstracting from matrices and thinking of L, R, S as (noncommutative) algebraic objects, 
their elimination from the above system leads to the mcKP hierarchy with product modified 
by Q (°f- 0S])- To some extent the above Riccati system thus expresses the mcKP hierarchy as 
a hierarchy of ordinary differential equations. 

Finite-size matrix Riccati equations, in particular with constant coefficient matrices as above, 
were discussed in a context related to integrable systems already long ago (see e.g. (49J 150]). 
but apparently not in the context of the KP hierarchy. A special infinite-size matrix Riccati 
system involving the shift operator in infinite dimensions appeared, however, in the framework 
of the Sato theory (see e.g. \51\ I13j). In the one-component case with B = I, the above Riccati 
system, with suitable conditions imposed on the coefficient matrices, also generates solutions of 
the BKP and the CKP hierarchy |52| . The Riccati system indeed generates solutions of various 
integrable systems and therefore deserves to be studied in its own right. 
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Remark 3. For fixed r G N, r > 1, and for some fixed B G B, let us consider the condition 
(HB) r Zo = Z P, 

with an N x N matrix P (over A). For the solution (|5.3p of the linear matrix system (|5.2p . this 
implies (HB) nr Z = ZP n for n G N, hence B nr (R nr X +Q nr Y) = XP n and B nr (S nr X+L nr Y) = 
YP n , and thus the algebraic Riccati equations B nr (S nr + L nr <$>) = YP n X~ 1 = <^XP n X^ 1 = 
QB nr (R nr + Q nr &). The corresponding equations (|6.2p of the Riccati system then read 

$t B , nr = B nr (S nr + L nr $) - $B nr (R nr - Q nr $) = 0, n = 1, 2, . . . . 

Hence <^> solves the (r, B)-reduction, i.e. the r-reduction (multicomponent version of rth Gelfand- 
Dickey hierarchy) with respect to B. 

7 Conclusions 

Any solution of a multicomponent Burgers (mcBurgers) hierarchy is a solution of the corre- 
sponding multicomponent KP (mcKP) hierarchy. Furthermore, there is a functional equation 
that determines the mcKP hierarchy and has the form of an inhomogeneous mcBurgers hierar- 
chy functional equation. We have also shown that the mcKP linear system is equivalent to 
a mcBurgers hierarchy, where the dependent variable has the structure of an infinite Frobenius 
companion matrix (which in particular makes contact with [9]). 

Moreover, we have shown how solutions of a mcKP hierarchy are obtained from solutions 
of a multicomponent linear heat hierarchy via a generalized Cole-Hopf transformation. An 
important subcase generates solutions of a mcKP hierarchy from solutions of a matrix linear 
system and we presented some explicit solution formulae. They comprise in particular Davey- 
Stewartson dromions and lump solutions. 

There is certainly a lot more to be (re)discovered using the rather simple but quite general 
method in Section [5] to construct exact solutions, but we are far from a systematic way to 
explore the properties of solutions obtained in this way. Furthermore, we have stressed the role 
of a matrix Riccati hierarchy in this context. 
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